A consistent statistical treatment of the renormalized mean - field t-J model 
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A variational treatment of the Gutzwiller - renormalized t-J Hamiltonian combined with the 
mean-field (MF) approximation is proposed, with a simultaneous inclusion of additional consistency 
conditions. Those conditions guarantee that the averages calculated variationally coincide with 
those calculated from the self-consistent equations. This is not ensured a priori because the effective 
Hamiltonian contains renormalization factors which depend explicitly on the mean-field averages. 
A comparison with previous mean-field treatments is made for both superconducting (d-RVB) and 
normal states and encompasses calculations of both the superconducting gap and the renormalized 
hopping amplitudes, as well as the electronic structure. The C^-symmetry breaking in the normal 
phase - the Pomeranchuk instability (PI) - is also analyzed. 

PACS numbers: 71.27.+a, 74.72.-h, 71.10.Fd 



t-J mode!= is regarded to reflect some of the essential 
physics of strongly correlated copper 3d x 2_ y 2 states in 
high-temperature superconductors^ In this model, the 
correlated hopping of electrons reduces strongly their 
band energy, so the latter, for the doping x < 0.1, be- 
comes comparable to the real-space pairing part induced 
by the kinetic exchanged However, the analytical solu- 
tions of t-J model are limited to the special cases for the 
one-dimensional system^ Under these circumstances, 
we have to resort to either exact diagonalization, 5 which 
is limited to small cluster systems or to the approximate 
methods. The latter include renormalization group, 6 
variational approach based on the Gutzwiller - projected 
wavefunctions (either treated within Monte Carlo tech- 
niques or by Gutzwiller approximation^) and various 
versions of the slave - boson approach^ Each of these 
methods seizes some of the principal features of these 
quasi- two-dimensional correlated states, although no 
coherent picture has emerged as yet. 

In this paper we concentrate on the Gutzwiller renor- 
malized mean-field (MF) theory for the t-J Hamilto- 
nian and formulate a variational procedure, with the ad- 
ditional conditions ensuring the self-consistency of the 
whole approach. Implementing such procedure is essen- 
tial (if not indispensable) for obtaining reliable results of 
the MF type. It is reassuring that some of the quantities 
such as the RVB gap magnitude or the hopping correla- 
tions (bond-parameter) do not change appreciably with 
respect to the earlier results j& whereas the others, such 
as the single-particle electronic structure, are altered re- 
markably. Furthermore, we illustrate the basic nontriv- 
iality of our approach on the example of the so-called 
Pomeranchuk instability discussed recently^ 

We start with the t-J model in its simplest form^ 



operators and the model parameters have the standard 
meaning.— 

To proceed further, effective mean-field renormalized 
Hamiltonian is introduce d 7 » 8 ' 10 ~ 13 i 15 i 16 which is taken in 
the following form 

H = X (*W'5*j4 CT C i<7 +H.C.) -fl^C^Cia 

(ij)a i<? 

~ l J a9ij(Xjich c io + Hc - - \Xij\ 2 ) 

- ^|j^(A ij -4c|_ (7 + H.c.-|A^| 2 ). (2) 

In the above expression, c! (cjo-) are ordinary fermion 
creation (annihilation) operators, Xij = ( c L c io-)' an< ^ 
Aij = {d-oCjtr) = (cj- a Ci a ) are respectively, the hop- 
ping amplitude (bond-parameter) and the RVB gap pa- 
rameter, both taken for nearest neighbors (ij). The 
renormalization factors gjj and gfj result from the 
Gutzwiller ansatz. The exchange part (Si ■ Sj) has been 
decoupled in the Hartree-Fock-type approximation and 
incorporates as nonzero all above bilinear averages ob- 
tained according to the prescription 



K C> 7 -> A s A t + A s A t - A s A t , (3) 
where t — t(K, 7) etc. and for any operator A 



A = (A)=Tr[Ap], 



(4) 



n, 



Pi 



with p being the density matrix for the mean-field Hamil- 
tonian to be determined. By taking the step from ([1} to 
([2]) we introduce essentially a non-Hartree-Fock-type of 
approximation, which differs from ([3]) due to the presence 



„t 



(ij) 



Sj—fi\^n ia ^P, (1) °f 9ij an d 9ij factors. Therefore, we may not be able use 

i rr p ff flip rlpnQitv cri&riAiTir nf tfip fnrm n = 7,~ ^ p~ @H 



where P labels the Gutzwiller projector which guarantees 
that no doubly occupied sites are present. The projected 



e.g. the density operator of the form p = Z 1 e~ 
Z = Tr[e~^ H ], as a proper grand-canonical trial state in 
the frame of variational principle based on the Bogoli- 
ubov inequality^ since then the self-consistency of the 



2 



approach (expressed by Eq.Q) may be violated. This is 
the reason, why in most of the previous mean held treat- 



ments, e.g 



801-13 



, the standard procedure encompasses 



diagonalizing of the bilinear Hamiltonian @, and subse- 
quently solving of the self-consistent (s-c) Bogoliubov-de 
Gennes (BdG) equations for Xij, Ajj an d In effect, 
this procedure does not refer to any variational scheme. 

The solution based solely on the s-c BdG equations, 
although acceptable, may not be fully satisfactory. This 
is because in the present situation we build up the entire 
description on the basis of MF Hamiltonian and hence we 
should proceed in a direct analogy to the exact (non-MF) 
case. Namely, our approach is based on the maximum 
entropy principle.— Such starting point provides us with 
a general variational principle, which may differ from that 
of Bogoliubov and FeynmanJ^ In other words, the value 
of the appropriate functional is minimized, with the self- 
consistency of the whole approach being preserved at the 
same time.— 

To tackle the situation, we define an effective Hamilto- 
nian H\ containing additional constraints, that is of the 
form 



H X = H 



i <y 

- E (^(i^-Xii)+H.c) 

- E (A£(c 4 _ CT c> - A tj ) + H.c. 



(5) 



where the Lagrange multipliers A-"', A*-, and A^ play 
the role of molecular fields. Moreover, the parameters 
Xij, A,;i, and rii = y~L(cJ„CinA coincide with those which 



appear in the renormalization factors <?*- and gf-, and 



which are taken in the form 



8.13 



iXjXjjl - Xj)(l - Xj) 



(1 - xf)(l - xj) + 8(1 - x lXj )\ Xl3 \ 2 + 16\Xij\ 4 



4(1 - Xi )(l- Xj ) 



(1 - xj )(1 - xj) + SxiXj&p) + 16/3+ (4) ' 



(6) 



with Xi = l- n,-, Pf 3 {n) = |A^|™ ± \ Xij \ n . 

When solving the model on a square lattice and in the 
spatially homogeneous case, there appear thus five mean 
fields, A ee (n,Xx,Xy, A y ), with Xr = Xtj,^ T = 
V2Ay, {{ij)\\T, t = x,y; as well as the same num- 
ber of the corresponding Lagrange multipliers, A ee 
(A, A*, A*, A*, A£), where A? = A*- , A^ = V2A£. Both A 

and A are assumed to be real. Apart from that, for given 
n we have to determine the chemical potential (i. The 
first step is the diagonalization of H\ via Bogoliubov- 
Valatin transformation, which yields 

Hx = E £k(7£o7ko + 7^7ki) + - E *) + C > ( ? ) 



with £ k = x/Sl + Dl, D k = V2E r £, rCOs(fc r ), and 
a = -2 J2 T T T cos(fc T ) - fi - A. Also, 

T T = -t lT g\ T + - a J t9 j tX t + A*, D T = \j T g J T A T + A* 

(8) 

X = An+ E (§-W(2;£+A?)+4x r A*+2A T A£). (9) 



For the sake of simplicity, we have included only the hop- 
ping between the nearest neighbors, although the gener- 
alization to the case with more distant hopping does not 
pose any principal difficulty. We define next the general- 
ized Landau functional, J- = — fi^ 1 ln(Tr[e - ^ x ]), which 
here takes the form 

F{A, A) = C + E ((& - E k ) - - In (l + e"^*)) , (10) 

k ' 

with inverse temperature f3 = 1/ksT. The equilibrium 
values oi A — Aq, X = Xq are the solution of the set of 
equations 



Va^ = 0, V a J" = 0, 



(11) 



for which (fTU)) reaches its minimum. This step is 
equivalent to the maximization of the entropy with 
the constraints^ Also, the grand potential ft and the 
free energy F are defined respectively as Q(T, V, fi) — 
F(T,V,MA (T,V,fj,),X (T,V,n)), and F = O + M iV. 
Note, that by taking the derivatives with respect to A 
only, and subsequently putting A = 0, the results reduce 
to the standard BdG self-consistent equations. 

Even though the present method can be regarded as 
natural within the context of statistical mechanics, to 
the best of our knowledge, it has not been utilized, in 
the form presented here, in the context of condensed 
matter physics problems. Also, in this respect, our ap- 
proach unifies individual features of the self-consistent 
variational MF treatments developed earlier— ~—, which 
in the T = limit can be obtained as particular cases. 
Parenthetically, the present method, together with the 
Gutzwiller approximation, provides also a natural justi- 
fication of some aspects of the slave-boson saddle-point 
approach, as some of the constraints coincide in both 
methods. 

We solve numerically first the system of equations (fTTj) 
on the lattice of A — 128 x 128 sites, using the pe- 
riodic boundary conditions and taking the parameters 
Jx = Jy = J — 1) tx = t y = — 3J, and for low tempera- 
ture k B T/J = 0.002 for the filling n = 7/8 = 0.875. Both 
the d-wave superconducting resonating valence bond (d- 
RVB) and the isotropic normal (N) solutions are ana- 
lyzed. The self-consistent variational results (denoted as 
var) obtained here and those obtained from BdG equa- 
tions are compared in Tables I and II. One sees that 
our value of the low-temperature free energy (per site), 
(c.f. Table I) in the d-RVB phase is slightly better than 



3 



the previous estimates albeit not much (—1.3661 for 
var, as compared to —1.3647 for s-c). It is slightly 
higher than that of the Variational Monte Carlo, which is 
Evmc/J = —1-3671, c.f. 13 Also the isotropic staggered- 
fiux (SF) phase has been found unstable against N state 
within both methods at this filling. In Table II we 
display microscopic quantities characterizing each solu- 
tion in that case and compare them with those obtained 
within standard s-c treatment. The differences are more 
pronounced for the RVB state. 

Table I. Comparison of the values of the 
thermodynamic potentials (per site). fl (F) stands for 
fl — XN (ft + fiN) for var and fl s - c (fl s -c + fJ, s - c N) for 
s-c methods, respectively. 



Therm. Pot. 


var RVB 


s-c RVB 


var N 


s-c N 


ft/A 


-5.75856 




-6.25862 




ft/A 


-1.07648 


-1.03575 


-1.11823 


-1.08025 


F/A 


-1.36614 


-1.36471 


-1.2955672 


-1.2955671 



Table II. Values of chemical potentials and MF 
parameters, fl stands for A + /x (var), and for /i s _ c (s-c) 

Variable var RVB s-c RVB var N s-c N 



A 

Xx = Xy 

A* 

£ z= = 

V2 y^L 



_Ay_ 

A 



5.01989 
-5.35094 
-0.33105 
0.18807 
-0.16985 
0.13199 

-0.01111 



5.67206 

-5.87473 
-0.37595 -0.20267 -0.24608 
0.19074 0.20097 0.20097 

-0.18369 
0.12344 0.00000 0.00000 

0.00000 
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FIG. 1: (Color online) Dispersion relations along the main 
symmetry lines in the Brillouin zone for a square lattice, of 
the size A x = A y — 128, and for the filling n — 0.875. Left (a): 
d-RVB solutions, right (b): N solutions. Triangles - earlier 
self-consistent results, circles - the present method. 

For the parameters listed in Tables I and II we have com- 
puted the quasiparticle energies in both the d-RVB and 
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the N states. Those are shown in Fig. la-b. The solid cir- 
cles represent our results, whereas the previous ones^ are 
drawn as triangles. The energy-dispersion reduction in 
our case is connected with presence of the constraints and 
results in a decrease of the bandwidth, which, in turn, is 
regarded as a sign of enhanced electron correlations. 

After testing the feasibility of our approach for fixed 
doping x, we now discuss systematic changes appearing 
as the function of shown in Figs. [5] and [3] 
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FIG. 2: (Color online) Doping dependences of the free energy 
(F V ar, F s -c), the chemical potentials /i, /i s _ c as well as that 
of A and A + fl for the d-RVB state, both within the present 
(var) and the standard (s-c) methods. 
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FIG. 3: (Color online) Left: 
bond-order parameters \x = 



0.1 0.2 0.3 
x = 1-n 



Doping dependence of the 
Xy, the superconducting or- 



der parameters A a 



-A H , their renormalized counterparts 



g t %x = g t Xy an( l <7 t A :E = —g t A y , as well as (right) of the 
quantities 2T X = 2T y and \/2D x = -y/2D v of Eq. ©, both 
for the s-c (triangles) and the var (circles) methods. 

We emphasize, the chemical potential /i is the first 
derivative of F/A with respect to n (c.f. Figj2]), unlike in 
some of the previous mean-fields treatments&ii~— (c.f. 
however Refi 10 i 14 ). This is also the reason why we dif- 
ferentiate between \x and fi = ji + A, even in the case of 
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the spatially homogeneous solution. The doping depen- 
dence of other relevant MF quantities is shown in FigJ3] 
The results are again close to those obtained from the 
BdG procedure, except for T r> (FigJ3Jb)), which enter 
the quasiparticle energies. 

So far we have focused on MF solutions with the sym- 
metry between x and y directions on the square latticed 
However, a spontaneous breakdown of this equivalence of 
the x- and y- directed correlations is possible already in 
the normal phase and is called the Pomeranchuk instabil- 
ity (PT)r- that manifests itself by lowering of the discrete 
Ci V symmetry of the Fermi surface. 
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FIG. 4: (Color online) Doping dependence of bond-order pa- 
rameters Xx and Xy (left), and the free-energy differences AF 
(right) both for x-y symmetric (N) and the Pomeranchuk (x- 
y symmetry-broken) states (PI) within both the present (var, 
filled circles) and the standard (s-c, triangles) methods, re- 
spectively. The vertical line marks the phase transition within 
the s-c method. For details, see main text. 

In Fig. |4] (a) the doping dependence of the bond- 
order parameters Xx an d Xy are displayed for the x- 
y symmetric (N) and the symmetry-broken (PI) solu- 



tions, both within our {x% ar ) and the standard (Xr _ °) 
methods. Within the s-c scheme, PI solution is found 
up to x w 0.091. However, a comparison of the re- 
spective free-energy differences, AF S _ C = F^_ c — F^J C 
and AF var = F v N ar - F™, (cf. Fig. g] (a)) reveals 
that this solution becomes unstable against N state for 
x « 0.021, thus the phase transition is certainly discon- 
tinuous. On the other hand, within our variational treat- 



ment the PI solution does not exist for x > xZ 



0.044, 



where AF var w 0, in qualitative agreement with what 
is expected for the continuous phase transition. From this 
analysis it is clear that the two methods of approach (s-c, 
var) yield qualitatively different predictions for PI. 

In summary, we have introduced self-consistency con- 
straints required within the variational mean-field ap- 
proach to the Gutzwillcr-renormalized mean-field t-J 
model. Such consistency conditions are indispensable 
from the basic statistical-mechanical point of view. Un- 
dertaking such a step results in consistent evaluations 
of the thermodynamic quantities, which in the present 
method are determined from the generalized Landau 
functional. A detailed comparison with the standard 
mean-field solution based on Bogoliubov-de Gennes self- 
consistent equations (i.e. that without constraints) is 
provided. Our method introduces quantitative and, in 
some cases, even qualitative corrections to the standard 
mean-field results. Other mean-field states such as flux 
phases or antiferromagnetism can be treated in the same 
manner. 
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